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Preface 


The  purpose  of  this  study  was  to  model  the  effect  charged  particle  radiation 
has  on  the  working  parameters  of  metal-oxide-semiconductor  field-effect  transis¬ 
tors  (MOSFETs).  There  are  two  individuals  I  would  especially  like  to  thank.  One 
is  my  thesis  advisor  Michael  Sabochick,  for  the  helpful  advice  and  guidance  he 
provided  throughout  the  project.  The  other  is  Roger  Tallon  of  the  Air  Force 
Weapons  Laboratory,  who  originally  suggested  this  as  a  possible  thesis  topic,  and 
provided  an  (at  the  time)  unpublished  copy  of  their  experimental  results. 
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Abstract 


A  simple  model  to  characterize  the  effects  of  heavy  charged  particles  on 
metal-oxide  semiconductor  field-effect  transistors  (MOSFET)  is  presented.  The 
model  is  applied  to  experimental  results  provided  by  the  Air  Force  Weapons 
Laboratory,  and  an  attempt  made  to  simulate  saturation  phenomena  observed  in 
the  threshold  voltage  change.  The  model  assumes  all  trapped  holes  are  within  a 
few  nanometers  of  the  oxide-semiconductor  interface,  and  takes  into  account  the 
resultant  counter  electric  field,  and  its  effect  on  the  yield  fraction  escaping  recom¬ 
bination.  An  equation  relating  threshold  voltage  change  as  a  function  of  dose  is 
derived  and  used,  but  does  not  duplicate  the  saturation  effect.  This  is  because 
charge  trapped  near  the  interface  reduces  the  internal  oxide  electric  field  very  lit¬ 
tle.  However,  similarities  in  p-channel  and  n-channel  MOSFET  damage  sensitivi¬ 
ties  indicate  considerable  bulk  charge  trapping.  This  would  modulate  the  internal 
electric  field  considerably.  To  determine  whether  this  accounts  for  saturation  will 
necessitate  refinement  of  the  model,  taking  hole  transport  and  bulk  trapping  into 


account. 


MODELING  THE  EFFECTS  OF  HEAVY 


CHARGED  PARTICLES  ON  METAL-OXIDE- 
SEMICONDUCTOR  FIELD  EFFECT  TRANSISTORS 


I.  Introduction  and  Background 


One  of  the  most  important  devices  used  in  very- large-scale  integrated  circuits 
is  the  metal-oxide-semiconductor  field-effect  transistor  (MOSFET).  Figure  1 
shows  a  schematic  of  an  n-channel  MOSFET.  The  device  works  as  follows. 

When  a  positive  voltage  is  applied  to  the  gate,  a  strong  electric  field  sets  up 
through  the  gate  oxide  and  penetrates  into  the  p-type  Si  substrate.  In  the  region 
just  beneath  the  gate  oxide,  this  field  repels  the  majority  carriers  (which,  for  n- 
channel  MOSFETs,  are  holes)  and  attracts  minority  carriers  (electrons),  thus 
forming  an  inversion  layer.  With  a  potential  difference  applied  between  source 
and  drain,  current  flows  through  this  inversion  layer  and  the  device  is  considered 
"on".  The  operation  of  a  p-channel  MOSFET  is  similar,  only  the  p-type  and  n- 
type  materials  are  reversed,  and  a  negative  voltage  is  applied  to  the  gate  contact. 

MOSFETs  are  very  sensitive  to  the  effects  of  ionizing  radiation.  The 
sequence  of  events  associated  with  irradiation  of  an  n-channel  MOSFET  is  shown 
in  Figure  2.  Ionizing  radiation,  such  as  gamma-rays,  electrons,  and  heavy  charged 
particles,  initially  create  many  electron-hole  pairs  throughout  the  gate  oxide.  A 
certain  fraction  of  these  pairs  are  then  annihilated  through  recombination.  The 
positive  applied  voltage  rapidly  sweeps  the  remaining  electrons  from  the  oxide  and 
into  the  metal  gate  contact.  The  holes  transport  more  slowly  toward  the  gate- 
semiconductor  interface,  where  a  certain  fraction  are  caught  in  deep  hole  traps 
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Figure  1.  Schematic  of  an  N-channel  MOSFET  (Ref  1:22) 


and  the  rest  disappear  into  the  semiconductor.  These  trapped  holes  permanently 
alter  the  electric  parameters  of  the  device,  primarily  by  lowering  the  threshold 
gate  voltage  required  to  create  the  inversion  layer.  In  fact,  if  enough  holes  are 
trapped  at  the  interface,  an  inversion  layer  may  form  without  applying  any  vol¬ 
tage  at  all.  For  p-channel  MOSFETs,  the  process  is  similar,  except  the  holes 
migrate  toward  the  metal-gate  contact  and  are  trapped  at  that  interface. 

Recent  experiments  have  been  conducted  by  the  Air  Force  Weapons  Labora¬ 
tory  to  characterize  the  effects  of  ionizing  radiation  on  MOSFETs.  References  3 
and  4  describe  this  work  in  detail.  Large  numbers  of  non-radiation-hardened 
Intersil  3N161  p-channel  and  3N171  n-channel  discrete  MOSFETs  were  irradiated 
with  protons,  electrons,  and  Co-60  gamma  rays.  After  exposure  to  specified  levels 
of  radiation  (under  various  bias  conditions),  the  transistor  gate  threshold  voltages 
were  measured  "in  situ",  and  the  changes  in  threshold  voltage  found  as  a  function 
of  radiation  angle  of  incidence  (the  angle  between  the  proton  track  and  the  electric 
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Figure  2.  Sequence  of  Events  Associated  with  Irradiation  of  MOSFET  (Ref 
2:67) 


field),  applied  field,  particle  energy,  and  total  ionizing  dose. 

Of  particular  interest  were  the  results  of  the  proton  irradiations.  In  this  case, 
the  MOSFETs  were  irradiated  at  room  temperature  with  protons  of  energies  from 
2-16  MeV.  The  angle  between  the  proton  track  and  the  electric  field  was  varied, 
and  data  collected  for  the  cases  of  0,  -45,  and  80  degrees.  Figures  3  and  4  show 
the  damage  sensitivity  (change  in  threshold  voltage  per  dose)  for  the  p-channel 
and  n-channel  MOSFETs,  as  a  function  of  particle  energy  and  angle  of  incidence. 
In  general,  the  damage  increases  with  both  particle  energy  and  angle  of  incidence. 
This  is  as  expected,  as  will  be  explained  in  the  theory  section  of  this  paper. 

What  is  not  expected,  however,  is  the  anomaly  where  the  damage  caused  by 
the  45  degree  protons  is  greater  than  that  of  the  80  degree  protons  for  particle 
energies  less  than  6  MeV.  When  first  observed,  the  researchers  thought  that  this 
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Figure  3.  P-channel  MOSFET  Damage  Sensitivity  as  a  Function  of  Particle 
Energy  and  Angle  of  Incidence  (Ref  4:1209) 

anomaly  was  due  to  the  additional  material  (silicon  dioxide  passivation  and 
aluminum  metallization)  the  80  degree  protons  had  to  pass  through  before  reach¬ 
ing  the  gate  oxide.  If  these  additional  material  lengths  were  greater  than  the 
mean  proton  ranges,  then  the  number  of  80  degree  protons  reaching  the  gate  oxide 
would  be  reduced.  However,  the  actual  distance  traveled  by  the  protons  was 
found  to  be  6  to  7  times  less  than  the  mean  ranges  of  2-6  MeV  protons.  This 
implied  that  all  of  the  protons  reached  and  passed  through  the  gate  oxide,  regard¬ 
less  of  the  angle  of  incidence. 

In  a  further  attempt  to  understand  this  anomaly,  the  data  was  replotted  in  a 
different  format,  shown  in  Figures  5,  6,  and  7.  Figure  5  shows  the  change  in 
threshold  voltage  as  a  function  of  dose  for  the  various  angles  of  incidence.  For 
this  case,  the  gate  bias  potential  was  -5  Volts.  Figures  6  and  7  show  essentially 
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Figure  4.  N-channel  MOSFET  Damage  Sensitivity  as  a  Function  of  Particle 
Energy  and  Angle  of  Incidence  (Ref  4:1209) 


the  same  thing,  though  with  different  gate  bias  potentials.  On  the  basis  of  these 
graphs,  the  researchers  deduced  that  a  saturation  effect  was  occuring  within  the 
test  sample,  that  had  the  same  effect  as  reducing  the  internal  electric  field  during 
irradiation. 

To  explain  this  phenomena,  the  researchers  proposed  that  the  buildup  of 
trapped  holes  at  the  interface  resulted  in  a  counter  electric  field.  The  yield  frac¬ 
tion  of  holes  escaping  initial  recombination  increases  with  both  the  angle  of 
incidence  and  the  magnitude  of  the  electric  field.  Initially,  the  higher  yield  frac¬ 
tion  of  the  80  degree  protons  results  in  a  larger  charge  buildup  at  the  interface. 
However,  this  larger  buildup  offsets  the  electric  field  within  the  oxide,  decreasing 
the  fraction  of  holes  escaping  recombination.  Thus,  the  threshold  voltage  change 
due  to  80  degree  protons  saturates  before  the  45  degree  case. 
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Figure  5.  Change  in  Threshold  Voltage  as  a  Function  of  Dose  and  Angle  of  Ir 
cidence  for  Applied  Gate  Voltage  of  -5  Volts  (Ref  4:1212) 

The  goal  of  this  thesis  project  was  to  give  this  rough  qualitative  model  a 
more  precise  quantitative  formulation,  and  verify  that  the  counterfield  effect  actu¬ 
ally  explains  the  anomaly.  The  next  section  describes  the  stages  shown  in  Figure 
2  in  more  detail.  Assumptions  used  in  the  model  are  stated,  and  mathematical 
formula  presented.  The  following  section  presents  the  major  results  and  analysis, 
and  is,  in  turn,  followed  by  the  conclusion. 
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Figure  6. 


Change  in  Threshold  Voltage  as  a  Function  of  Dose  for  80  Degree 
Angle  of  Incidence  and  Gate  Voltages  of  -5  and  -10  Volts  (Ref  -4:1213) 
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II.  Theory 


This  section  describes  the  assumptions  and  mathematical  equations  that 
enter  into  each  stage  of  the  model.  Considered  in  turn  are  columnar  phenomena 
including  initial  electron-hole  formation  and  recombination,  bulk  phenomena 
including  electron  and  hole  transport,  and  interface  phenomena  including  deep- 
hole  trapping,  and  radiation-induced  traps.  Following  this,  an  equation  relating 
the  change  in  threshold  voltage  to  the  total  ionizing  dose  is  developed.  Both  the 
dependence  of  the  electric  field  on  trapped  charge,  and  the  dependence  of  yield 
fraction  escaping  recombination  on  electric  field  are  considered  in  this  develop¬ 
ment. 


Columnar  Phenomena 

As  the  protons  pass  through  the  oxide,  they  lose  energy  through  ionization 
and  excitation  of  electrons  across  the  energy  bandgap;  electron-hole  pairs  are 
formed.  The  amount  of  energy  needed  to  create  an  electron-hole  pair  in  silicon 
dioxide  has  been  estimated  to  be  about  18  eV  (Ref  5:1520).  Protons,  as  well  as 
other  heavy  charged  particles,  leave  a  very  dense  column  of  electron-hole  pairs  in 
their  track.  However,  to  determine  the  initial  charge  produced  by  an  incremental 
dose,  the  electron-hole  pairs  are  assumed  to  be  evenly  distributed  throughout  the 
bulk  of  the  oxide.  This  incremental  charge  is  given  as 

AC=KgADAlfv  (1 

where 


AC  =  charge  (Coulomb) 


Kg  =  generation  constant  (Coulomb/cm3— rad) 

AD  =  dose  (rad— Si02) 

A  =  area  of  oxide  (cm2) 

1  =  length  of  oxide  (cm) 
fy  =  fraction  yield 

The  generation  constant,  Kg,  is  found  from 

yr  _  (.OlJ/kg— rad)(1.609  X  10~19C/hole)p  ^ 

^  (1.609  X  1(T19)(1000  g/kg)W 

where 

p  =  density  (g/cm3) 

W  =  energy  per  electron— hole  pair  formed  (eV/hole) 

For  Si02  the  density  is  2.2  g/cm3  and  the  energy  required  to  form  an  electron  hole 
pair  is  18  eV.  Filling  these  numbers  into  Eq  (2)  gives  a  generation  constant  of 
1.222  X  10~6  C oulomb/cm3— rad. 

Immediately  after  the  electron-hole  pairs  are  produced  in  the  dense  column,  a 
certain  fraction  of  them  will  recombine.  The  fraction  that  survive,  fv.  is  a  func¬ 
tion  of  the  particle  type  and  energy,  the  electric  field  within  the  oxide,  and  the 
incident  angle  between  the  particle  track  and  electric  field.  To  find  the  fraction 
yield,  a  numerical  procedure  developed  by  Oldham  (described  in  references  6  and 
7)  is  used.  A  description  of  this  procedure  follows. 

Figure  8  shows  a  schematic  of  the  track  of  a  charged  particle  after  passing 
through  the  silicon  dioxide  insulator.  Looking  down  the  direction  of  the  track, 
the  initial  charge  density  configuration  is  assumed  to  be 


Figure  8.  Charged  Particle  Track  Through  the  Oxide 


n±  =  electron— hole  density  (cm-3) 

N0  =  initial  linear  density  (cm-1) 
b  =  Gaussian  radius  (cm) 
r  =  Distance  center  of  the  track  (cm) 

The  electrons  and  holes  then  move  according  to  the  equation  (Ref  7:2695) 

dn±  dn± 

-fa-  =  D±V"n±  +  M±EX-^—  -  Qn+n_ 

where 

D±  =  Diffusion  coefficient  (cur's) 

=  mobility  (cm2,V— s) 

Ex  =  Component  of  electric  field  in  the  x— direction  (V/cm) 
a  =  Recombination  coefficient  (cm3/s) 

The  first  term  on  the  right  hand  of  Eq  (-1)  represents  diffusion,  the  second  term 
represents  drift  under  the  influence  of  the  external  field,  and  the  third  term 
represents  recombination. 
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The  actual  program  (found  in  the  Appendix)  sets  up  a  grid  in  the  x  and  in¬ 
direction  with  z  in  the  direction  of  the  proton  track,  and  y  perpendicular  to  the 
electric  field  for  all  proton  angles  of  incidence.  The  initial  electron-hole  density  is 
found  from  Eq  (3)  for  each  point  in  the  grid.  These  densities  are  then  placed  into 
a  finite  difference  form  of  the  right  hand  side  of  Eq  (4),  and  dn±/dt  found  at  each 
point  on  the  grid.  The  new  densities  are  then  found  from 

5n . 

M^new)  =  n±(r-told)  +  ^j~At  (5) 


where  At  is  the  increment  in  time.  These  new  densities  are  inserted  into  the  right 
hand  side  of  Eq  (4),  and  the  process  repeated  until  the  electron  and  hole  columns 
separate  under  the  influence  of  the  external  field.  Integrating  n+(r,t)  over  the  area 
of  the  grid  gives  a  new  linear  hole  density  N(t). 

Figure  9  shows  the  electron-hole  columns  moving  past  each  other  under  the 
influence  of  an  electric  field  of  arbitrary  direction.  For  the  case  of  a  particle  track 
perpendicular  to  the  electric  field,  the  fraction  yield,  fy,  will  be  given  simply  by 
N(t)/N0,  where  N(t)  is  the  linear  density  after  the  columns  have  separated.  The 
analysis  up  to  this  point  is  only  good  for  the  region  between  the  dashed  lines  in 
Figure  9,  but  it  can  be  generalized  to  three  dimensions.  If  it  is  assumed  that  the 
hole  column  to  the  right  of  the  right-hand  dashed  line  remains  unchanged,  then 
the  fraction  yield  can  be  found  from 


f„  = 


T 

/  N(t)dt 


NCT 


(6) 


where  T  is  the  time  required  for  the  right  hand  of  the  electron  column  to  cross  the 
left  hand  of  the  hole  column. 


The  parameters  fj,±,  D±,  a,  N0,  and  b  in  Eqs  (3)  and  (4)  are  found  in  the  fol¬ 
lowing  way.  The  room  temperature  mobilities  of  silicon  dioxide  are  found  quite 


Figure  9.  Electron  and  Hole  Columns  Moving  Under  the  Influence  of  Arbitrary 
Electric  Field  (Ref  7:2696) 


readily  from  the  literature.  The  values  used  in  this  analysis  were 

(J+=2X  IO-5cm2/ V— sec  (Ref  8:2012)  and  -  20  cm2/Y-sec  (Ref  9:1333).  The 

diffusion  coefficients,  D±,  are  then  found  from  the  mobilities  by  using  the  Einstein 


relation  (Ref  7:2696) 


D±- 


H±kT 


where  k  is  the  Boltzmann  constant.  T  is  the  temperature,  and  e  is  the  electron 
charge.  The  recombination  coefficient.  a,  is  found  by  using  an  expression  derive 
by  Langevin  (Ref  7:2696) 

(/■*+  +  /Oe 


where  c  is  the  dielectric  constant  of  silicon  dioxide,  and  c0  is  the  permittivity  of 
free  space.  The  initial  linear  density,  N0,  depends  on  particle  energy  and  can  be 
found  from  stopping  power  tables  (Ref  10:310).  The  value  of  the  Gaussian  radius, 
b,  was  taken  from  Reference  7  where  it  was  used  as  a  variable  pir-inieter  to  fit 


this  recombination  model  to  experimental  data.  The  value  was  found  to  be  3.5 
nanometers  (Ref  7:2697). 

It  should  be  pointed  out  that  this  model  is  only  valid  for  situations  where  the 
electron-hole  pair  separation  distance  is  short  enough  that  the  pairs  lose  their 
identity.  In  this  case  electrons  readily  recombine  with  holes  from  different  pairs. 
This  is  generally  true  of  the  particle  tracks  left  by  heavy  charged  particles, 
because  of  their  high  linear  energy  transfer.  For  electrons  and  gamma-rays,  the 
pair  separation  distance  is  large  enough  that  recombination  is  pretty  much  res¬ 
tricted  to  electron- holes  of  the  same  pair.  The  model  developed  to  handle  these 
cases  is  called  the  geminate  model  (Ref  1:27). 

In  general,  the  yield  fraction  fy  depends  on  three  factors.  First,  it  depends  on 
the  energy  of  the  particle.  This  is  because  the  initial  linear  density,  N0  in  Eq  (3). 
generally  decreases  with  increasing  particle  energy,  at  least  in  the  ranges  of 
interest  here.  The  lower  the  initial  density,  the  less  probability  there  is  that 
recombination  will  occur,  thus  the  yield  fraction  increases.  Second,  yield  fraction 
also  depends  on  the  strength  of  the  electric  field.  The  stronger  the  field,  the  more 
quickly  the  electron  and  hole  columns  separate,  decreasing  the  time  available  for 
recombination.  Thus,  yield  fraction  increases  with  increasing  electric  field.  Third, 
the  yield  fraction  depends  on  the  angle  of  incidence.  This  can  be  understood  by 
considering  the  two  extremes.  When  the  particle  track  is  perpendicular  to  the 
field,  the  two  columns  are  separated  rapidly,  which  increases  yield  fraction.  When 
the  particle  track  is  parallel  to  the  field,  the  columns  are  pulled  through  each 
other,  enhancing  recombination,  and  reducing  the  yield  fraction.  Thus,  yield  frac¬ 
tion  increases  with  increasing  angle  of  incidence. 

This  increase  in  yield  fraction  with  particle  energy,  electric  field  and  angle  of 
incidence  is  the  cause  of  the  general  characteristics  of  the  curves  seen  in  Figures  3 
through  7.  In  these  figures,  the  change  in  threshold  voltage  also  increases  with 


particle  energy,  electric  field,  and  (except  for  particle  energies  below  6  MeV)  angle 
of  incidence.  As  will  be  seen  below,  the  change  in  threshold  voltage  is  directly 
proportional  to  the  charge  accumulated  at  the  interface  of  the  oxide  and  thus  the 
yield  fraction. 

Bulk  Phenomena 

After  escaping  initial  recombination  the  remaining  electrons  are  rapidly  swept 
from  the  oxide  by  the  applied  gate  field.  This  process  takes  a  few  picoseconds  to 
complete.  Two  assumptions  are  made  about  this  stage.  First,  no  electrons  are 
trapped  within  the  oxide.  Second,  having  escaped  columnar  recombination,  the 
electrons  do  not  recombine  further  with  holes  previously  trapped  in  the  oxide. 

.After  the  electrons  are  swept  from  the  device,  the  holes  transport  toward  the 
interface  much  more  gradually.  .Assume  all  the  holes  move  with  the  same  velo¬ 
city,  v  =  fj,+ E.  For  an  oxide  of  length  280  nanometers  and  an  electric  field  of 
1  X  103  Y/cm  at  room  temperature,  it  takes  a  hole  1.40  X  10~5  seconds  to  cross 
the  oxide.  This  time  is  short  enough  to  justify  ignoring  any  transient  effects 
caused  by  transport.  To  simplify  the  analysis  then,  all  the  holes  produced  by  an 
incremental  dose,  AD  (that  escape  recombination),  are  assumed  to  instantaneously 
transport  to  the  interface.  Bulk  trapping  within  the  oxide  is  ignored. 

Interface  Phenomena 

.As  the  holes  pass  through  the  interface  region,  a  certain  fraction,  fj,  are 
caught  in  deep  hole  traps  within  about  10  nanometers  of  the  interface.  The  frac¬ 
tion  trapped  depends  on  the  temperature  and  electric  field,  but  is  mostly  a  result 
of  the  processing  techniques  used  during  manufacture.  The  interface  region  has 


high  local  strain  and  a  deficiency  of  oxygen  atoms.  This  gives  rise  to  a  large 
number  of  Si-Si  bonds  instead  of  the  usual  Si-O-Si  bonds  found  in  silicon  dioxide. 
When  a  hole  passes  one  of  these  bonds,  it  can  break  the  bond  and  recombine  with 
one  of  the  bonding  electrons.  This  results  in  a  Si  atom  retaining  the  remaining 
electron  from  the  broken  bond,  and  the  positive  charge  residing  with  the  other  Si 
atom  (Ref  1:34). 

Since  the  trapping  fraction  depends  primarily  on  the  manufacturing  process, 
it  can  vary  quite  radically  between  different  MOSFETs.  Generally,  fT  can  be  as 
little  as  1-2  percent  for  hard  oxides,  10-20  percent  for  good  quality  commercial 
oxides,  and  as  much  as  50-70  percent  in  very  soft  commercial  oxides  (Ref  1:33). 

The  spatial  dependence  of  the  charge  trapped  at  the  interface  is  also  depen¬ 
dent  on  the  relative  hardness  of  the  device.  Hard  oxides  have  distributions  that 
fall  off  exponentially  as  one  moves  away  from  the  interface  into  the  oxide.  On  the 
other  hand,  the  softer  oxides  have  distributions  that  remain  relatively  constant 
(Ref  11:1207). 

These  trapped  holes  can  persist  from  a  few  hours  to  years.  Annealing  does 
take  place,  in  which  electrons  from  the  Si  semiconductor  tunnel  into  the  oxide  and 
recombine  with  the  trapped  holes.  However,  since  the  time  frame  for  annealing  is 
considerably  longer  than  the  experiment,  it  is  ignored  in  this  model. 

This  model  assumes  that  a  certain  fraction,  fT,  of  the  holes  escaping  recombi¬ 
nation  will  be  trapped  in  a  constant  charge  distribution,  p ,  extending  a  distance  X 
from  the  interface.  The  incremental  charge  density  due  to  an  incremental  dose  is 
then 


■v = 


(9) 


in  the  region  near  the  interface  and  0  elsewhere.  X  is  presumably  about  10  nm 
(Ref  1:33).  The  trapping  fraction  is  uncertain,  but  since  these  devices  are  reported 


as  non-radiation  hardened  (Ref  3:4393),  it  is  probably  greater  than  10  percent. 
This  quantity  is  left  undetermined  in  the  analysis,  but  is  assumed  constant. 

Radiation-induced  interface  traps  exist  very  close  to  the  interface,  and  within 
the  Si  band  gap.  Their  occupancy  is  determined  by  the  position  of  the  Fermi 
level  at  the  interface,  and  can  be  negative  (usually  with  n-channel  devices),  posi¬ 
tive  (usually  with  p-channel  devices),  or  neutral.  In  either  case,  the  presence  of 
these  trapped  charges  results  in  shifts  to  the  threshold  voltage.  At  the  present 
time,  the  precise  mechanisms  responsible  for  the  formation  of  these  states  are  not 
well  understood,  and  a  subject  of  considerable  debate  (Ref  1:37).  These  states  are 
difficult  to  model.  Since  the  deep  hole  traps  are  dominant  (especially  at  the  lower 
dose  levels  applicable  here),  interface  traps  are  ignored  in  this  analysis. 


Threshold  Voltage  vs  Dose 

In  this  section,  an  analytical  expression  relating  the  change  in  threshold  vol¬ 
tage  to  the  total  ionizing  dose  is  found,  so  that  the  curves  in  Figures  5  through  7 
can  be  simulated.  This  expression  will  include  the  dependence  yield  fraction  has 
on  electric  field  and  incident  angle. 

The  change  in  threshold  voltage  for  arbitrary  charge  density  is  given  by  (Ref 
12:199) 

AVfh  =  J^Mdx  (10) 

o  f 

where 

1  =  the  length  of  the  oxide  (crn) 
p(x)  =  density  of  holes  (cm-3) 
e  =  permittivity  of  Si02 

Assuming  the  charge  density  is  constant  between  I  —  X  and  1.  and  zero  elsewhere. 
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the  threshold  voltage  change  becomes 


A  _  Ptx2  -  2M) 

th  2e 


(H) 


or 


d  =  2X1) 

tn  2e 


(12) 


where  dAVth  is  the  incremental  change  in  threshold  voltage  resulting  from  an 
incremental  change  in  charge  density.  Substituting  the  value  for  A p  given  by  Eq 
(9)  into  Eq  (10)  gives 


dAVlh  «  K&fyfTl  (X-  21  ^  AD  (13) 

In  Eq  (13),  everything  is  assumed  constant  except  the  yield  fraction.  fv,  which 
depends  on  the  electric  field,  which  in  turn  depends  on  the  charge  distribution. 

We  next  derive  an  expression  relating  the  internal  electric  field  to  p  and  thus 
AVth.  Figure  10  shows  a  schematic  cross  section  of  a  MOSFET,  with  an  areal 
density  Qx<  located  at  arbitrary  x'  within  the  oxide.  This  charge  induces  a  charge 
density  Qo  at  the  gate,  and  Qj  at  the  semiconductor.  Ignoring  any  difference  in 
work  function  between  the  metal  gate  and  semiconductor,  it  can  be  shown  that 


Qo  = 


(U) 


and 

Qi  =  ~^Qx 


(13) 


Using  Gauss’  Law,  the  internal  electric  field  due  to  charge  density  Qx.  becomes 

E(x)  =  --!1— tp,  (16’ 

for  0  <  x  <  x'  and 


Semiconductor 


Figure  10.  Schematic  Cross  Section  of  MOSFET  With  Arbitrary  Areal  Charge 
Density  Qx,  at  x' 


E(x)  =  4-Qx, 


for  x'  <  x  < 


Eqs  (16)  and  (17)  can  be  generalized  to  an  arbitrary  charge  distribution  p(x) 
as  shown  in  Figure  11.  Taking  into  account  the  field  due  to  the  external  gate  vol¬ 
tage  and  letting  Qx  =  p(x'  )dx' ,  the  electric  field  at  x  becomes 

e(x)  =  /Artodx-  _ M  +  \  (1S, 


Expanding  the  second  term  on  the  right  and  rearranging  gives 


T 


Figure  11.  Schematic  Cross  Section  of  MOSFET  with  Arbitrary  Charge  Densi 
tv 


The  first  term  on  the  right  hand  side  of  Eq  (19)  is  merely  —  AVth/l.  Concern 
ing  ourselves  with  the  field  only  at  x  <  1  —  X.  the  second  term  becomes  —pX/t. 


Finding  p  in  terms  of  AVth  from  Eq  (11)  and  inserting  into  this  second  t 


erm  gives 


E(x)  = 


*th  X 
1  21  -  X 


For  simplicity,  the  yield  fraction  is  assumed  linear  with  electric  field  and  can 


be  given  the  form 


fv(E)  =  a,E  +  b, 


where  a*  and  b^  will  depend  on  the  angle  of  incidence  for  a  proton  of  given  energy 
Substituting  Eq  (20)  into  Eq  (21)  will  give  fv  as  a  function  of  change  in  threshold 
voltage.  Placing  fy(AVth)  into  Eq  (14)  yields  an  equation  of  the  form 


+  .*  V  /  ,*  ,»  ,■  .• 


,  •  -  - .  • .  v  v  ' j  ■r-  • '  s'  v  -/ ■: 


III.  Results  and  Analysis 


The  previous  section  developed  the  model  used  to  simulate  the  effect  charged 
particle  radiation  has  on  MOSFET  devices.  This  section  describes  the  results  of 
applying  this  model  to  the  actual  experiment  described  in  the  introduction  and 
background  section. 

Figure  12  shows  the  yield  fraction  as  a  function  of  particle  energy  for  angles 
of  incidence  of  0,  45,  and  80  degrees  using  the  columnar  recombination  model.  An 
electric  field  of  1.79  X  105  V/cm  was  used.  This  field  was  found  by  dividing  the 
gate  voltage  (5  Volts)  by  the  thickness  of  the  oxide  (280  nanometers).  The  fraction 
yield  increases  with  both  particle  energy  and  angle  of  incidence,  and  the  charac¬ 
teristics  of  these  curves  are  similar  to  those  shown  in  Figures  3  and  4  except  for 
the  data  anomaly. 

In  finding  a  functional  relationship  between  the  change  in  threshold  voltage 
and  the  dose  (accounting  for  oxide  field  modulation  due  to  hole  trapping  near  at 
the  interface)  it  was  assumed  that  the  fraction  yield  was  a  linear  function  of  the 
electric  field  within  the  oxide.  Figure  13  shows  the  fraction  yield  as  a  function  of 
electric  field  for  4  MeV  protons  with  angles  of  incidence  of  45  and  80  degrees. 

This  figure  was  found  by  using  the  columnar  recombination  model.  As  can  be 
seen,  the  fraction  yield  increases  with  both  electric  field  and  angle  of  incidence, 
but  the  curves  are  not  linear  throughout  the  range  of  interest.  To  properly  use 
Eq  (23).  a  linear  approximation  was  taken  from  0.1  MY/cm  to  0.25  MV/cm  (the 
general  range  of  the  oxide  electric  field  during  the  irradiation).  The  coefficients  in 
Eq  (21)  are  found  to  be  a#  =  4.11  X  10~7  cm/V  and  b#  =  .0539  for  an  angle  of  45 
degrees,  and  =  4.88  X  10  7  cm/\  and  b$  =  .0650  for  an  angle  of  80  degrees. 

I  sing  these  values  in  Eq  (23)  with  X  =  10  nm,  1  =  280  nm.  and  V  =  5  Volts, 
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Figure  12.  Fraction  Yield  as  a  Function  of  Particle  Energy  and  Angle  of  In¬ 
cidence. 

AVth  was  found  as  a  function  of  dose  for  trapping  fractions  of  .20.  and  .10.  The 
results  are  shown  in  Figure  13.  .As  expected,  the  threshold  voltage  rises  with  dose, 
angle  of  incidence,  and  trapping  fraction.  However,  the  sought  after  saturation 
effect  is  missing. 

The  reason  the  saturation  effect  is  absent  is  straightforward.  If  all  of  the 
charge  is  assumed  to  accumulate  within  10  nm  of  the  interface,  the  charge  giving 
rise  to  a  threshold  voltage  change  equivalent  to  the  gate  voltage  will  (from  Eq 
(20))  only  decrease  the  internal  field  to  0.982  of  its  original  value.  This  is  gen¬ 
erally  true  for  charges  at  the  interface— they  modulate  the  internal  electric  field 
very  little,  certainly  not  enough  to  cause  tlm  *=•> t u ration  effect  seen  in  Figures  '> 
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Figure  13.  Fraction  Yield  as  a  Function  of  Electric  Field  and  Angle  of  Incidence 

through  7. 

This  raises  the  question  whether  all  of  the  charge  is  actually  trapped  near  the 
interface.  There  is  evidence  that  it  is  not.  Figures  3  and  4  show  the  damage  sen¬ 
sitivity  for  both  n-channel  and  p-channel  MOSFETs.  Under  applied  bias,  holes 
created  by  radiation  accumulate  at  the  gate  interface  in  p-channel  MOSFETs.  and 
at  the  semiconductor  interface  in  n-channel  MOSFETs.  From  Figures  3  and  4.  it 
can  be  seen  that  the  damage  sensitivity  is  a  little  bit  greater  in  the  n-channel  dev¬ 
ices.  but  not  much.  Our  original  assumption  was  that  the  charge  was  trapped 
within  lOnm  of  the  interface.  Using  Eq  (10).  a  given  charge,  p,  at  the  interface  of 
the  p-channel  device  results  in  a  threshold  voltage  change  of  —30 p/e.  If  this  same 
charge  is  distributed  at  the  semiconductor  interface  of  the  n-channel  device.  Eq 
(10)  gives  a  AVth  =  —  2750p/e,  which  is  55  times  greater.  Of  course  this  assumes 


Figure  14.  Threshold  Voltage  Change  As  a  Function  of  Dose.  Angle  of  Incidence, 
and  Trapping  Fraction 


the  trapping  fractions  are  equal  at  both  interfaces,  which  is  unlikely.  However,  it 
seems  improbable  that  the  trapping  fractions  are  that  different.  If  a  charge  p  is 
evenly  distributed  throughout  the  bulk  of  the  oxide,  then  the  change  in  threshold 
voltage  is  equivalent  for  both  types  of  MOSFETs. 

Thus  it  appears  that  there  is  a  significant  amount  of  charge  trapped  in  the 
bulk  of  the  oxide.  The  hole  transport  process  is  actually  very  dispersive.  Some 
holes  move  out  of  the  oxide  quickly,  while  others  take  a  long  time  to  do  so.  The 
assumptions  made  about  hole  transport  trapping  are  generally  true  for  thin  oxides 
less  than  1 00  nm  and  high  electric  fields,  greater  than  1  MV/cm.  For  the  oxides 
considered  here,  the  thickness  was  280  nm  and  the  field  1.79  X  105  \  /cm.  With 
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these  parameters,  it  is  likely  that  a  significant  proportion  of  charge  is  trapped  in 
the  bulk  of  the  oxide  (Ref  13:3940). 

The  implication  of  large  bulk  trapping  can  be  seen  from  Figure  15.  This 
figure  shows  the  electric  field  as  a  function  of  distance  for  a  charge  density  distri¬ 
buted  evenly  throughout  the  oxide.  The  charge  density  is  assumed  large  enough 
to  give  rise  to  a  threshold  voltage  change  equal  to  the  applied  gate  voltage.  From 
Eqs  (11)  and  (19)  the  internal  electric  field  is  given  as 

E(x)-2Vg-L  (24) 

so  E(0)  =  0  and  E(l)  =  2V_/I.  The  distributed  charge  causes  greater  modulation 
throughout  the  bulk  of  the  oxide,  than  that  charge  accumulated  at  the  surface. 

To  the  left  of  the  midpoint  in  Figure  15  the  electric  field  falls  significantly, 
causing  a  reduction  in  yield  fraction.  On  the  other  hand,  the  field  rises 
significantly  to  the  right  of  the  midpoint,  which  increases  the  yield  fraction. 

From  Figure  13,  the  yield  fraction  seems  to  slope  off  with  decreasing  electric  field, 
so  the  decline  in  fraction  yield  in  the  left  hand  region  may  dominate  the  increase 
in  the  right  hand  side,  resulting  in  an  overall  reduction  in  fraction  yield.  Whether 
this  reduction  is  great  enough  to  cause  the  saturation  effect  remains  questionable. 

It  could  well  be  a  contributor,  but  there  are  other  things  to  consider. 

For  instance,  it  was  assumed  earlier  that  the  electrons,  having  escaped  colum¬ 
nar  recombination,  were  swept  from  the  oxide  without  further  recombination. 

F’or  low  hole  densities  and  high  electric  fields  this  assumption  is  certainly  valid. 
However,  the  cross  section  for  electron-hole  recombination  is  proportional  to  E-'  2 
(Ref  14:3203).  The  reduced  field  in  Figure  15  could  result  in  significant  recombi¬ 
nation  in  the  region  to  the  left  of  midpoint.  This  could  he  another  source  of  the 
saturation  effect,  especially  if  the  field  reaches  zero  within  the  oxide.  Then  the 
electrons  would  stop  and  recombine  with  a  hole  (Ref  16:1522-1523). 


Figure  15.  Electric  Field  Within  the  Oxide  from  an  Evenly  Distributed  Charge 
Density 

A  third  assumption  made  in  this  mode!  was  that  the  fraction  of  holes 
trapped  at  the  interface  was  constant.  However,  it  is  known  that  the  hole  trap¬ 
ping  cross  section  is  dependent  on  the  electric  field  and  varies  with  E-1  *  (Ref 
15:1192).  The  increasing  field  at  the  interface,  as  seen  in  Figure  15  may  result  in 
a  decline  in  trapping  fraction,  that  could  also  be  a  contributor  to  the  saturation 
effect. 

These  three  field  modulated  mechanisms— yield  fraction,  electron  recombina¬ 
tion,  and  hole  trapping  may  all  contribute  to  the  saturation  effects  seen  in  Figures 
5-7.  To  verify  this  will  require  a  more  refined  model  that  not  only  takes  these 
mechanisms  into  account,  but  also  accounts  for  the  wide  lispersion  inherent  in  the 
hole  transport.  This  will  be  required  to  characterize  the  hole  buildup  in  the  bulk 


IV.  Conclusion 


In  conclusion,  this  simple  model  has  shown  that  the  electric  field  reduction 
due  to  holes  trapped  at  the  interface  cannot  cause  the  saturation  effect.  Charges 
trapped  close  to  the  interface  cannot  reduce  the  internal  electric  field  to  any  great 
degree— certainly  not  to  the  extent  required  to  effect  the  yield  fraction  of  holes 
escaping  columnar  recombination.  However,  similar  p-channel  and  n-channel  dam¬ 
age  sensitivities  indicate  that  considerable  charge  is  trapped  within  the  bulk  of  the 
oxide.  Charge  distributed  within  the  oxide  bulk  modulates  the  electric  field  con¬ 
siderably.  This  modulated  electric  field  could  contribute  to  the  saturation  effect  in 
three  ways—  reduction  of  the  yield  fraction  throughout  the  oxide,  increase  in  elec¬ 
tron  recombination,  reduction  of  hole  trapping  at  the  interface. 

To  verify  this,  the  present  model  would  have  to  be  refined  in  a  number  of 
ways.  First,  the  dispersive  nature  of  hole  transport  would  have  to  be  considered 
so  that  reasonable  estimates  of  the  trapped  hole  distribution  (both  at  the  interface 
and  within  the  bulk  of  the  oxide)  can  be  made.  Second,  Eqs  (10)  and  (19)  would 
have  to  be  modified  to  numerically  find  the  threshold  voltage  and  electric  field  for 
an  arbitrary  charge  density.  Third,  the  nonlinearity  between  yield  fraction  and 
electric  field  would  have  to  be  accounted  for,  so  that  overall  yield  fraction  due  to 
a  widely  varying  internal  electric  field  may  be  found.  Fourth,  the  electron  tran¬ 
sport  and  electric  field  dependent  recombination  cross  section  would  have  to  be 
introduced.  Finally,  the  hole  capture  cross  section  dependence  on  electric  field 
would  have  to  be  taken  into  account.  Introducing  these  refinements  into  the 
model  should  allow  determination  of  the  primary  cause  of  the  saturation  effect 
and  explain  the  data  anomaly. 


Appendix:  Columnar  Recombination  Computer  Code 


This  appendix  contains  the  code  used  to  find  the  yield  fraction  of  electrons 
and  holes  escaping  inita.il  recombination.  The  following  algorithm  is  used: 


Step  1. 


Step  2. 


Step  3. 


Step  4. 


Step  5. 

Step  6. 
Step  7. 
Step  8. 
Step  9. 
Step  10. 
Step  11. 

Step  12. 
Step  13. 


Input  the  initial  linear  density,  the  mobilities,  the  diffusion  coefficients, 
the  recombination  coefficient,  the  electric  field,  the  angle  of  inci 
incidence,  the  column  radius,  and  the  length  of  the  oxide  (lines  22-23). 

Set  up  the  x  and  y  grid.  In  this  case  the  maximum  values  in  the  x  and 
y  directions  and  the  number  of  grid  spaces  are  read  in.  The  program  is 
set  up  to  take  advantage  or  the  symmetry  in  the  y  direction  (lines  30- 
45). 


Input  and  initialize  time  related  elements.  These  include  the  time  incre¬ 
ment  and  the  ratio  of  time  increments  between  each  step.  This  ratio  is 
inserted  to  help  shorten  the  run  times.  The  time  variable  is  initialized 
and  the  maximum  time  calculated  (lines  51-53). 

Set  up  boundary  conditions  on  the  grid.  In  this  routine  the  boundaries 
are  assumed  large  enough  that  the  electrons  and  holes  separate  well 
before  the  columns  reach  the  boundary  time  so  the  electron  and  hole 
densities  at  xmax,  xmin,  and  ymax  are  all  set  to  zero  (lines  58-73). 

Determine  the  linear  density  of  particle  at  time  equal  to  zero  (lines  108- 
124). 


Determine  the  yield  fraction  for  each  time  step  (lines  132-216). 
Determine  new  time  and  time  ratio  (lines  133-134). 

Set  new  boundary  condition  along  ymin  (lines  139-1  12) 
Determine  the  new  density  at  each  point  on  grid  (lines  144-183). 
Find  the  diffusion  term  in  Eq  (2)  (lines  151-160). 


Find  the  drift  term  in  Fq  (2)  (lines  164-170). 

Find  the  recombination  term  in  Eq  (2)  (line  17  1). 

find  the  time  rate  of  change  of  the  holes  electrons  (lines  179- 
180). 


Step  14. 


f  ind  the  new  charge  density  at  the  new  time  (lines  188-196). 


Step  15.  Find  the  number  of  charge  carriers  per  unit  length  and  thus  the 

new  yield  (lines  200-212). 

Step  16.  Print  out  the  time, electron  and  hole  density,  and  the  yield  (line 

213). 

Appendix  of  variables. 

alpha  -  recombination  coefficient  in  em3/seo 
b  -  Gaussian  radius  in  nanometers 
bcm  -  Gaussian  radius  in  centimeters 

check  -  variable  used  to  determine  whether  the  point  is  outside  the  column  or  not 
delx  -  the  length  of  each  element  in  the  grid  in  the  x-  direction  in  nanometers 
dely  -  the  length  of  each  element  in  the  grid  in  the  y-  direction  in  nanometers 
delt  -  the  time  increment  in  seconds 

derele(i.j)  -  time  rate  of  change  of  electron  charge  density  in  /cm3*s 
derhol(i.j)  -  time  rate  of  change  of  hole  charge  density  in  /cm3*s 
diffn  -  diffusion  term  in  the  electron  density  rate  of  change  equation 
diffp  -  diffusion  term  in  the  hole  density  rate  of  change  equation 
dndx  -  first  derivative  with  respect  to  x  of  electron  density 
dpdx  -  first  derivative  with  respect  to  x  of  hole  den  density 
dplus  -  diffusion  coefficient  of  holes  in  cm2/s 

drift n  -  drift  term  in  the  electron  density  rate  of  change  equation 
driftp  -  drift  term  in  the  hole  density  rate  of  change  equation 
dneg  -  the  diffusion  coefficient  of  the  electrons  in  cm2/s 
d2ndx2  -  second  derivative  with  respect  to  x  of  the  electron  density 

d2ndy2  -  second  derivative  with  respect  to  y  of  the  electron  density 

d2pdx2  -  second  derivative  with  respect  to  x  of  the  electron  density 

d2pdy2  -  second  derivative  with  respect  to  y  of  the  hole  density 


V**.* 


I 


e  -  the  electric  fiefield  applied  across  the  oxide  in  Volts/cm 
eleden(i,j)  -  the  electron  charge  density  in  charges  /cm3 

eledenx  -  the  density  of  electrons  along  the  x  axis 

eledeny  -  the  density  of  electrons  along  the  y  axis 

elenum  -  the  number  of  electrons  along  the  y  axis  in  charge/cm 

elenuml  -  the  number  of  electrons  above  the  y  axis  in  charge/cm 
elenumo  -  the  total  number  of  electrons  equal  to  elenum+elenum  1  *2 
hx  -  thethe  length  of  each  element  in  cm  in  the  x  direction 
hy  -  the  length  of  each  element  in  cm  in  the  y  direction 
holden(i,j)  -  the  hole  charge  density  in  charges  /cm3 
holdenx  -  the  hole  density  along  the  x  axis  in  charges/cm3 
holdeny  -  holthe  hole  density  along  the  y  axis  in  charges/cm3 
holnum  -  the  number  of  holes  along  the  y  axis  in  charges/cm 
holnuml  -  number  of  holes  above  the  y  axis  in  charges/cm 
holnumo  -  total  number  of  holes  equal  to  holnum-f  holnum  1  *2 
length  -  length  of  oxide  in  cm 

m  -  number  of  elements  between  upper  and  lower  boundaries  in  the  grid  in  x 
direction 

midx  -  midpoint  of  the  grid  in  the  x  direction 
midy  -  midpoint  of  the  grid  in  the  y  direction 
muneg  -  mobility  of  electrons  in  cm2/V*sec 
muplus  -  mobility  of  holes  in  cm2/V*sec 

n  -  number  of  elements  between  the  upper  and  lower  boundaries  of  the  grid  in  the 
y  direction 

no  -  initial  line  density  in  charges/cm 

recom  -  recombination  term  in  the  density  rate  of  change  equation 

thetad  -  angle  of  incidence  between  the  particle  track  and  the  electric  field  in 
degrees 

thetar  -  angle  of  incidence  in  radians 
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time  -  time  for  current  step 

tmax  -  maximum  time  for  program  to  run 

tration  -  ratio  of  delt  from  one  iteration  to  that  of  the  previous  one 

x(i)  -  values  of  x  for  each  grid  point  i  in  nano-  meters 

xmax  -  maximum  value  of  x  in  nanometers 

xmin  -  minimum  value  of  x  in  nanometers 

y(i)  -  values  of  y  for  each  grid  point  i  in  nano-  meters 

yield  -  fraction  yield  for  columnar  recombination 

vmax  -  maximum  value  of  y  in  nanometers 

ymin  -  minimum  value  of  y  in  nanometers 


0001  c  This  code  is  designed  to  find  the  concentration  of  holes  and  elec- 

0002  c  trons(in  charges/cm**3)  and  the  yield  fraction  of  holes  escaping 

0003  c  recombination  aue  to  interaction  of  charged  particles  (alphas  and 

0004  c  protons)  with  silicon  dioxide. 

0005 

0006  implicit  none 

0007  real  pi,no,muneg,muplus,dneg,dplus.alpha.e,thetad,thetar 

0008  real  b,bcm,xmax,xmin,ymax,ymin,delx,dely,x(20l),y(102) 

0009  real  delt,eleden(201, 102), holden(201, 102), check. yield 

0010  real  derele(201,102),derhol(201,102),hx,hy,t  max,  time,  length 

0011  real  elenum.holnum 

0012  real  elenumo,holnumo,eledenx(201,100),holdenx(201.100) 

0013  real  eledeny(l01,100),holdeny(101.100),elenum  1  .holnum  1  ,d2ndx2 

0014  real  d2ndy2,diffn,dndx,driftn,recom,tratio 

0015  real  d2pdx2.d2pdy2,diffp,driftp,dpdx 

0016  integer  n,m,midx,midy,i,j 

0017  open(  unit=  1, file=' rec. dat', st  at  us= 'old '.READONLY) 

0018 

0019 

0020  c  This  section  initializes  a  number  of  the  constants. 

0021 

0022  read!  l,*)no,muneg,muplus,dneg, dpi  us, alpha 

0023  read(l,*)e,thetad,b, length 

0024  pi =3. 1 416 

0025  thetar=thetad*pi/l80 

0026 

0027 

0028  c  This  section  sets  up  the  x  and  v  grid. 

0029 

0030  read(l,*)xmax,xmin,ymax,vmin,m.n 

0031  bcm=b*lE-07 

0032  midx=m/2-(-l 

0033  midv=2 

0034  delx=(xmax-xminh/m 

0035  delv=(ymax-ymin)/n 

0036  hx=delx*  IE-07 

0037  hy=dely*  IE-07 

0038  x(l)=xmin 

0039  y(l)=ymin 

0040  do  10  i=2.m+l 

0041  x(i)=x(i-l)+delx 

0042  10  continue 

0043  do  20  j=2,n+l 

0044  y(j  )= v(j- 1  )-|-del  v 

0045  20  continue 

0046 

0047 

0048  c  This  section  inputs  the  time  related  elements  such  as  the  time 

0049  c  increment,  and  the  maximum  time  considered. 

0050 

0051  read(l,*)delt,tratio 

0052  tmax=(iength/cos(thetar))/(muneg*cos(thetar)*e) 

0053  time=0.0 

0054 


0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 


c  This  section  inputs  the  boundary  conditions  for  the  charge  density 
c  at  xmin  and  xmax. 

do  30  j=2,n+l 
eledenfl,j)=0.0 
eleden(n+l,j)=0.0 
holden(l,j)=0.0 
holden(n+ 1  ,j  )=0.0 
30  continue 

c  This  section  inputs  the  boundary  conditions  for  the  charge  density 
c  at  ymax. 

do  40  i=2,m 

eleden(i,n+l)=0.0 
holden(l,n+l)=0.0 
40  continue 


c  This  section  initializes  the  charge  densities  at  time  equal  zero. 

do  50  i=2,m 
do  60  j— 2,n 

check=-(x(i)**2+y(j)**2)/b**2 
if(check.lt.-100.00)tnen 
eleden(i,j)=0.0 
holden(i  ,j)=0.0 
else 

eleden(i,j)=(no/(pi*bcm**2))*exp(-(x(i)**2+y(j)**2)/b**2) 

holden(i,j)=eleden(i,j) 

endif 

60  continue 
50  continue 


c  This  section  finds  the  charge  density  along  the  x  and  y  axis  at 
c  time  equal  to  zero. 

do  70  i  — l,m  +  l 

eledenx(i.l)=eleden(i,midy) 

holdenx(i,l)=holden(i,midy) 

70  continue 

do  80  j  =  1  ,n  +  l 

e!edeny(j,l)=eleden(midx.j) 

holdeny(j,l)=holden(midx.j) 

80  continue 


c  This  section  determines  the  total  number  of  charge  carriers(per 
c  unit  length)  at  time  equal  to  zero. 


elenum=0.0 


0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 


holnum=0.0 
elenuml=0.0 
holnuml=0.0 
do  90  i=2,m 
do  100  j=3,n 

elenuml=elenuml+eleden(i,j)*(hx*hy) 

holnuml=holnuml+holden(i,j)*(hx*hy) 

100  continue 
90  continue 

do  110  i=2,m 

elenum=elenum+eleden(i,2)*hx*hy 
holnum=holnum-l-holden(i,2)*hx*hy 
1 10  continue 

elenumo=elenum-l-elenuml*2 
holnumo=holnum+holnuml*2 
yield— holnumo*  del  t 
print*,time,elenumo,holnumo,yield 
print*, 'max  time=’.tmax 

c  This  section  determines  the  time  derivative  of  the  charge  density 
c  the  new  charge  density  and  the  total  charge  as  a  function  of  time. 

120  if(time.gt.tmax)go  to  130 
time=time+delt 
delt=delt*tratio 
elenum=0.0 
holnum=0.0 
elenuml=0.0 
holnuml—  0.0 
do  140  i  =  l,m  +  l 

eleden(i,l)=eleden(i,3) 

holden(i,l)=eleden(i,3) 

140  continue 

do  150  i=2,m 
do  160  j=2,n 


c  This  section  finds  the  diffusion  term  for  the  electrons  from  finite 
c  difference  form. 

d2ndx2=(eleden(i+l,j  )-2*eleden(i,j)+eleden(i- 1  ,j  ))/hx*  *2 
d2ndy2=(eleden(i,j  +  I  )-2*eleden(i,j  )+eleden(i.j- 1  ))/hy**2 
diffn=dneg*(d2ndx2+d2ndy2) 

c  This  section  finds  the  diffusion  term  for  the  holes  from  finite 
c  difference  form. 

d2pdx2=(holden(i  +  l  ,j  l-2*holden(i,jl+holden(i- 1  ,j)l/hx** 
d2pdy2=(  holden(i,j  + 1  )-2*  holdenf i.j  )-*-liolden(  i.j- 1  ))/hy*  * 
diffp=dplus*(d2pdx2+d2pdy2) 

r  This  section  finds  the  drift  term  for  the  elect  r 
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0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 


dndx=(eleden(i  +  l,j)-eledenfi-l,j))/(hx*2) 
driftn=muneg*e*sin(thetar)*dndx 

This  section  finds  the  drift  term  for  the  holes. 

dpdx=(holden(i  +  l.j}-holden(i-l.j))/(hx*2) 
driftp  =  muplus*e*sin(thetar)*dpdx 

This  section  finds  the  recombination  term  for  both  holes  and  electrons. 

recom=-alpha*eleden(i,j)*holden(i,j) 

This  section  then  finds  the  time  rate  of  change  of  the  electron  and 
hole  density  for  each  x  and  y. 

derele(i,j)=diffn+driftn-f  recom 
derhol(i,j)=diffp+driftp+recom 

continue 

continue 

This  section  finds  the  new  charge  density  at  each  point  for 
new  time. 

do  170  i=2,m 
do  180  j=2.n 

eleden(i,j)=eleden(i.j)-Hlerele(i.j)*de)t 

holden(i.j)=holden(i.j)+derhol(i.j)*delt 

if(eleden(i.j).lt.0.0)then 

eleden(i,j)=0.0 

endif 

continue 

continue 

This  section  finds  the  number  of  charge  carriersfper  unit  length). 

do  190  i=2,m 
do  200  j=3.n 

elenum  1  =elenum  1  +eleden(i.j)*(hx*hy) 
hoi  num  1  =  holn  um  1  +  holden(i  .j  )*  ( hx  *  hy ) 
continue 
continue 
do  210  i=2.m 

elenum=elenum+eletlen(i,2)*hx*hy 
holnum  =  holnum  +  holden(i.2)*hx*hy 
continue 

elenum=elenuml*2+elenum 
holnum  =  holnuml  *2+holnum 
yield=yield+holnum*delt 
print*,  time.elenum.  hoi  num.  yield 

go  to  120 
continue 
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